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The equilibrium positions of the multi-particle classical Calogero-Sutherland-Moser 
(CSM) systems with rational/trigonometric potentials associated with the classical root sys¬ 
tems are described by the classical orthogonal polynomials; the Hermite, Laguerre and Jacobi 
polynomials. The eigenfunctions of the corresponding single-particle quantum CSM systems 
are also expressed in terms of the same orthogonal polynomials. We show that this interest¬ 
ing property is inherited by the Ruijsenaars-Schneider-van Diejen (RSvD) systems, which are 
integrable deformation of the CSM systems; the equilibrium positions of the multi-particle 
classical RSvD systems and the eigenfunctions of the corresponding single-particle quantum 
RSvD systems are described by the same orthogonal polynomials, the continuous Hahn (spe¬ 
cial case), Wilson and Askey-Wilson polynomials. They belong to the Askey-scheme of the 
basic hypergeometric orthogonal polynomials and are deformation of the Hermite, Laguerre 
and Jacobi polynomials, respectively. The Hamiltonians of these single-particle quantum 
mechanical systems have two remarkable properties, factorization and shape invariance. 


§1. Introduction 


Exactly solvable and/or quasi-exactly solvable multi-particle quantum mechan¬ 
ical systems have many remarkable properties. Especially, those of the Calogero- 
Sutherland-Moser (CSM) systems^^^^^ and their integrable deformation, called the 
Ruijsenaars-Schneider-van Diejen (RSvD) systems,"^)’have been extensively stud¬ 
ied. They have many attractive features at both classical and quantum mechanical 
levels and are related to topics studied in other fields. 

The equilibrium positions of the classical CSM systems with rational/trigono- 
metric potentials associated with the classical root systems are given by the zeros 
of the classical orthogonal polynomials, namely, the Hermite, Laguerre and Jacobi 
(Chebyshev, Legendre, Gegenbauer) polynomials.®^^®^ For example, the Hamilto¬ 
nian of the H-type multi-particle classical Calogero system is 
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This equation can be converted into a differential equation for f{y) = nj^=i(y “ Vj)- 

It is easy to show that ^ -^ const. Thus the polynomial f{y) satisfies 

f" — 2yf' = Af, in which the constant A is determined as A = —2n. This is the 
differential equation for the Hermite polynomial Hn{y). Therefore the equilibrium 
positions are given by the zeros of the Hermite polynomial. On the other hand, the 
corresponding single-particle quantum system is described by the Hamiltonian of 


H = -I- + const. 

2m 2 


( 3 ) 


This is just a harmonic oscillator, and its eigenfunctions are given by (ground state 
eigenfunction) x (Hermite polynomials). Interestingly, the equilibrium positions 
of the multi-particle classical system and the eigenfunctions of the corresponding 
single-particle quantum system are described by the same orthogonal polynomials, 
the Hermite polynomials. This is not a fortuitous coincidence. In fact, the H-type 
(or C-type) Calogero systems and the HC-type Sutherland systems have the same 
property. In these cases, the relevant orthogonal polynomials are the Laguerre and 
Jacobi polynomials, respectively. 

In this article, we establish that the above interesting property is inherited by 
the RSvD systems, which are integrable deformation of the CSM systems. Namely, 
we show that the equilibrium positions of the multi-particle classical RSvD systems 
with rational/trigonometric potentials associated with the classical root systems and 
the eigenfunctions of the corresponding single-particle quantum RSvD systems are 
described by the same orthogonal polynomials; the continuous Hahn (special case), 
Wilson and Askey-Wilson polynomials, which are deformation of the Hermite, La¬ 
guerre and Jacobi polynomials, respectively. These polynomials are members of the 
Askey-scheme of the basic hyper geometric orthogonal polynomials. Deforma¬ 

tion patterns are shown in Table I. For the trigonometric An-i systems, the situation 
is rather trivial; the equilibrium positions are equally-spaced and translationally in¬ 
variant, and the single-particle system is a free theory. Although we are able to relate 
them to Chebyshev polynomials of the first kind,®)’^^)’^^) we do not demonstrate 
this here. 


Table I. Deformation patterns. 


potential 

root system 

CSM system 


RSvD system 

rational 

A„-i 

Hermite poly. 


continuous Hahn poly. 

rational 

BC„ 

Laguerre poly. 


Wilson poly. 

trigonometric 

BC„ 

Jacobi poly. 


Askey-Wilson poly. 


For these single-particle quantum mechanical systems, the Hamiltonian has two 
remarkable properties, factorization and shape invariancef^^^^^^ Shape invariance 
is an important ingredient of many exactly solvable quantum mechanics. In our case, 
the shape invariance determines the eigenfunctions and spectrum from the data of 
the ground state eigenfunction and the energy of the first excited state. 

This article is organized as follows. In §2 we recapitulate the essence of the 
CSM systems with rational/trigonometric potentials associated with the classical 
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root systems. We briefly outline the argument that the equilibrium positions of the 
multi-particle classical systems and the eigenfunctions of the corresponding single¬ 
particle quantum systems are described by the same orthogonal polynomials; the 
Hermite, Laguerre and Jacobi polynomials. In §3 we generalize this argument to 
the RSvD systems. The Hermite, Laguerre and Jacobi polynomials are deformed to 
the continuous Hahn (special case), Wilson and Askey-Wilson polynomials, respec¬ 
tively. Section 4 is for a summary and comments. In order to study single-particle 
systems, we repeatedly use the shape invariance of the single-particle Hamiltonian. 
In appendix some useful techniques for shape invariant single-particle (‘discrete’) 
quantum mechanics are summarized along the idea of Crum.^®) 

§2. Calogero-Sutherland-Moser systems 

In this section we summarize the basics of the Calogero-Sutherland-Moser sys¬ 
tems with rational/trigonometric potentials associated with the classical root sys¬ 
tems. 

2.1. Models 

A CSM system is a multi-particle quantum (or classical) mechanical system gov¬ 
erned by a Hamiltonian H{p,q) (or classical one, The dynamical 

variables are the real-valued coordinates q = ^{qi, • • • , qn) and their canonically con¬ 
jugate momenta p = *{pi, -'' ^Pn)- For the quantum case we have pj = —ih-£p. We 
keep dimensionful parameters, e.g., mass, angular frequency, the Planck constant, 
etc. The coordinate qj has the dimension of length. 

The Hamiltonian of the CSM systems is 


Hcs {P, ^ ^P] + Fes (q), (4) 

where the potential Vcs{q) can be written in terms of the prepotential W{q), 



-|- h 


d'^W{q)\ 


( 5 ) 


The explicit forms of the potential Vcs(<?) and the prepotential W{q) are as follows: 
(i) rational : 


1 .0 " 


Vesiq) = ^ E 


9(9 - 1 ) 1 


1=1 


{qj -qkY 


— ^hwn{l -|- g{n — 1)) , 




W{q) =-^-mevq'^ + ^ 9h^og J ^ \qj - qk\, 


i=i 


l<j<k<n 


(6a) 

(6b) 
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(ii) rational BCn '■ 


n . 


i=i 


2„2 


-mta g • + 


{gs + 9L){gs + dL - 1) 


+ 




E 


f gjviigM - 1) , ghiigM - i) 


+ 


V (gj - gkY {qj + quY 

- hion{gs + 5L + 5 + gM{n - 1)) , 


(7a) 


n ^ 

W{q) = -Y^ -mujq] + ^ 


w\qj-qk\ + '^ogj^\qj + qk 


j = ^ 

n 


l<j<k<n 


+ ^^(ssftlog ^ IijI +9lfilog ^\2g,\) , 


i=i 

(iii) trigonometric BCn ■ 


Vcsiq) = 




TT 


2mL‘^ 
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{gs + gL){gs + g i - i) gdgi - i) 


i=i 

j-2 2 ^ 


sin2 f g,- 


cos2 iQj 


E 


gM{gM - 1) gM{gM - i) 


Vsin2 2 ^^^. _ g-j ^2 ^ 




2mL‘^ 


(^{gs + 2 ( 7 l + gM{n — 1 )) + g‘M\{n^ — 1 )^ , 

W{g)= E 5M^(log|sinf(gj-gfc)|+log|sinf(gj+gfc)|) 
l<j<k<n 
n 

+ '^{9sfi'^og\smj^qj \ + 5^/1 log|sin . 

i=i 


(7b) 


(8a) 


(8b) 


The constant terms in Vcsiq) 9-'^® the consequences of the expression @ in terms 
of the prepotential. A constant shift of Wiq) does not affect In the above 
formulas g,gs,gM and ql are dimensionless coupling constants and we assume they 
are positive. The other notation is conventional: m is the mass of the particles, 
uj is the angular frequency, h is the Planck constant (divided by 27r) and L is the 
circumference. All these parameters are positive. 

Note that the Hamiltonian can be expressed as a sum of factorized forms: 


j=i J J 


(9) 


Since the independent coupling constants are qm and gs + Ql, this BCn model is equivalent 
to the Bn model (or the C„ model). 
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2.2. Equilibrium positions of n-particle classical mechanics 

The classical Hamiltonian q) is obtained from the quantum one H{p, q) 

by the following procedure: (a) regard pj is a c-number; (b) after replacing the 
dimensionless coupling constants g, gi, 92 -,''' by the dimensionful ones g = gh,gi = 
gift, 92 = 52^) • ■ ■) assume 9 , 91 , 92 -,''' are independent of h ; (c) take the h ^ 0 
limit. In the same way, and are also obtained. 

The canonical equations of motion of the classical systems are 

dqj _ dH^^^^^{p,q) dpj _ dH^^^^^{p,q) 

dt 


dpj 


dt 


dqj 


The equilibrium positions are given by the stationary solution 

p = 0, q = q, 

in which q satisfies 

g^class(o^g) 


( 10 ) 

( 11 ) 


dqj 


= 0 . 


(j = l,...,n) 


( 12 ) 


q=q 


——— OV^class 

For the CSM system, m becomes — '^gq . \ = 0 and it is equivalent to the 


condition^^^ 


q=q 


g^class(^) 

dqj 


= 0 . 


(j = 1,... ,n) 


(13) 


q=q 


For the rational A^-i system, the above condition reduces to ©• Let us consider a 
polynomial whose zeros give the equilibrium positions: 


(i) rational A^-i : 

(ii) rational BCn '■ 

(hi) trigonometric BC„ 


/n(2/) = n(y-V^9T 

j=l V 

9M 

n 

fniO = • 




(14a) 

(14b) 

(14c) 




Then Ea. (ll3l) can be converted into a differential equation for fn{x), which is satisfied 
by the Hermite, Laguerre and Jacobi polynomials, respectively. The results are (see, 
for example. Refs. 6),8), 13), 22)) 

(i) : Uy)=H~{y), (15a) 

(ii) : /n(2/) = Li“)“(?/'), a = 


(hi): /n(0=4“’^^“(0 


a = 


9M 
9S + 9L 
9M 


- 1, (15b) 

- 1, /3 = |£ _ 1. (15c) 

9M 


(y), L(“)(y2)= ^ 


Here iL„(y) = 2"iL”°“= 
pioi,/3) momc are the Hermite, Laguerre and Jacobi polynomials, respectively. 
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2.3. Eigenfunctions of single-particle quantum mechanies 
The Hamiltonian 0 for single-particle case (n = 1 ) reads 


d? 1 2 2 

H = -H —muj X - TvjJ . 

2m dx^ 2 2 


(ii) : 

(iii) : 


d"^ 1 2 2 

H = —-+ -moj X 


H = - 


2m dx"^ 
2m dx^ 


, 5(5 - 1 ) fe / , 


+ 


2mL‘^ V sin^ j^x 


+ 


9 TT 

COS^ J^x 


fg{g-l) , gfg'-lh 


2mLf 


{ 9 + 9 


/\2 


(16a) 

(16b) 

(16c) 


where x = qi. These are a harmonic oscillator, a harmonic oscillator with a cen¬ 
trifugal barrier, and the Poschl-Teller potential,respectively, with some constant 
energy shifts. By introducing the dimensionless variable y, 


(i), (ii) : y = 


muj 


TT 


(iii) : 


these Hamiltonians are reduced to dimensionless ones 7i: 


(i) : 

H = hwH, 

n 

(ii) : 

H = hioH, 

n 

(iii) : 


n 


1 1 2 1 

'2d^ “^ 2 ^ ~ 2 ’ 

1 I I g(g-l) 

2dy2 “^2^ 2y2 


9 2 ’ 


1 , 9(9 - 1 ) 9'{9' - 1 ) _ / , / n 2'\ 

2 V dy'^ sin^ n cos2 v / 


sm y 


cos^ y 


(17) 

(18a) 

(18b) 

(18c) 


They have two properties, factorization (Hill) and shape invarianee (IA-5|) . In the 
notation of the appendix, the necessary data of the shape invariance are the in¬ 
variant (pre)potential W(y; A), the set of parameters A, the elementary shift of the 
parameters S and the energy of the first excited state Ti(A): 


(i) : W(y) = no A (no S), Si = 1, (19a) 

(ii) : W(y; A) =-h fflogy, A = 5 , <5 = 1, £:i(A) = 2 , (19b) 

(iii) : W(y ; A) = 5 log sin yg' log cos y , 

>^ = { 9 , 9 '), ^ = (1,1), SiW =‘^{9 + 9 '+ ^)- (19c) 


From the formulas (IAT6I) - ()AT7I) and (lATlII . we obtain the eigenfunctions Pn<Po and 
the corresponding eigenvalues Sn of 77: 


(i) : Pn{y) oc Hn{y), Sn = n, ( 20 a) 

(ii) : Pn{y,X)(xLn ^^(y^), ^^(A) = 2 n, ( 20 b) 

(iii) : Pniy;X)(xPn ^\cos2y), SniX) = 2n{n + g + g'). ( 20 c) 


Thus the equilibrium positions of the multi-particle classical systems and the 
eigenfunctions of the single-particle quantum systems are described by the same 
orthogonal polynomials; the Hermite, Laguerre and Jacobi polynomials. 
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§3. Ruijsenaars-Schneider-van Diejen systems 


In this section we generalize the results of the previous section to the Ruijsenaars- 
Schneider-van Diejen (RSvD) systems with the rational/trigonometric potentials as¬ 
sociated with the classical root systems. 

3.1. Models 

The RSvD system is an integrable deformation of the CSM system.The 
Hamiltonian of the RSvD systems is 


n 

H{p,q) = -mc^ ^ Vj {q) e^cPj ^ Vj (q) * + \JVj{qY e" Vj {q)-Vj {q)-Vj {qf^ , 

( 21 ) 


where Vj{q) is given by 


Vj{q) = w{qj) v{qj - qk) x 


k = l 


1 for An-i, 

v{qj + qk) for BCn . 


( 22 ) 


We use the conventional notation that Vj{q)* is the complex conjugate of Vj{q). 

. 1 zpi — 

Since the operators = e cause hnite shifts of the wavefunction in the 

imaginary direction f{q) = f{qi, • • • , qjpi:^, • • • , q-n)), we call these systems 

‘discrete’ dynamical systems. (Sometimes they are misleadingly called ‘relativistic’ 
version of the CSM. See Ref. 24) for comments on this point.) The basic potential 
functions v{x) and w{x) are as follows: 

(i) rational An-i ; 


v(x) = 1 — i 


■ k 9 


me X 


w{x) = (^ + i — + i — 3^ , 


(ii) rational BCn '• 


{ \ 1 -^90 

v[x) = 1 — i - 

me X 

. 


. a>2 


w{x) = (l-|-i — x][l + i — x)(l — i -) (1 — i 


. h 91 


me X 


h 


92 


mex-i 


h ! ’ 

2mc 


(hi) trigonometric BCn ■ 


v{x) = 


sm j^x 


(23a) 

(23b) 

(24a) 

(24b) 

(25a) 


w{x) = 


s[nf {x-i^gi) sinf(x-i^^^ 52 ) 

h 


sm 


cos iix- i—JY cos f (x - - i— 52 ) 


cos J^X 


cosf(x-i^ 


X 


(25b) 
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Here g, go, gi, 92, di and 92 are dimensionless coupling constants and c is the (ficti¬ 
tious) speed of light. We assume they are all positive. 

Note that the Hamiltonian m can be expressed as a sum of factorized forms, 

.. n _ _ _ 

H = -mc^ — \JVj{q)* ^Vj{q)* ^ . 

j=i 

(26) 

We also remark that in the c —> 00 limit, the RSvD systems reduce to the CSM 
systems, 

lim H{p,q) = Hcs{p,q), (27) 

c—^oo 

where the correspondence of the parameters is 

(i) : i0i+uj2 = i^, g = g, (28a) 

(ii) : uJi+u ;2 = u;, go = gn, 9 i + 92 = gs + 9 L , (28b) 

(hi): 90 = 9 M, 91 +92 =93 + 91 , g'l + 92 = 9 L ■ (28c) 


3.2. Equilibrium positions of n-particle classical mechanics 

The classical Hamiltonian H^^^^{p,q) (and also w^^^^{x)) is 

obtained by the prescription given in 112.21 For the RSvD system, the equation for 
the equilibrium positions m is equivalent to the condition^^) 

^^'"ss(g-) = Vf^^^iq)* > 0 . (j = 1,..., n) (29) 

This equation without the inequality sign can be rewritten in a Bethe-ansatz-like 
form, 


n 

k=i 


^class(^^. tjclassi'g^. + g^)* 


w 


class 




W 


class 


(Qj 


(j = 1 ,... ,n) 


(30) 


(For the An-i type systems, + q^) and -|- qk)* should be omitted.) 

Let us consider a polynomial whose zeros give the equilibrium positions: 


(i) rational An-i ■ 

fn{y)-fl{y J qj), 

j=l V 

(31a) 

(ii) rational BCn ■ 

n 

/„(!/)=n(!'^ "“‘«i). 

(31b) 

(iii) trigonometric BCn : 

n 

■ /n(0 = • 

i=i 

(31c) 

Then Ea.(l29l) can be converted into a 

functional equation for 

For the 


simplest case (i) with 102 = 0, this functional equation can be solved explicitly by 
using its generating function. However, such a straightforward solution of the func¬ 
tional equation cannot be easily obtained in the other cases. For this reason, we 
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solved numerically and guessed a recursion relation for {/„(x)}. We hypothe¬ 
size that the polynomials {fnix)} satisfy the three-term recurrence fn+iix) — {x — 
an)fn{x) + bnfn-iix) = 0 (n > 0, f-i{x) = 0, /o(x) = 1) for some constants 
and bn- The three-term recurrence implies that {fnix)} are orthogonal polynomials. 
Then we can show that the polynomials {fnix)} with this three-term recurrence 
solve the functional equation. The results are^^^’^^^ (see also Refs. 25)~27)) 


(i) : 

(ii) : 
(hi) : 


fniy) = P: 
fniy) = W, 


monic 

'n 


' / rrec^ mc^ 

mc^ 

rrec^ 

mc^ 

\ -Vi - •} 

. y uJig uJig 

CO 25 ’ 


’ xi2g 

9 9 

9 



, / me 9 me 

me 

9i 

92 \ 

1 — y 1 — •) 

Voiifiio <^i5o 

<^ 2 ff 0 

1 — ? 
9o 

9o^ ’ 

Zirsi 27r32 


0 -/ 

27rpj^ 


^ • g mcL ^ g mcL 

,-e 

mcL 

,-e 


27i-g2 


27rgQ . 

g mcL I . 


With proper 


normalization, these orthogonal polynomials are known as:^^^ 


(32a) 

(32b) 

(32c) 


(i) : the continuous Hahn polynomials, (33a) 

Pniy,ai,a2,bi,b2) = \in + ai + a2 + bi + b2- l)nPn°“''(2/; ai, 02 , ^ 2 ), 

re! 

(ii) : the Wilson polynomials, (33b) 

lTn(y^; ai, 02, as, 04) = (-l)”(re - l)n oi, 02,03,04), 

(hi) : the Askey-Wilson polynomials, (33c) 

Pnif.\ai,a2,a^,ai\q) = 2”(ai02a3a4g”"^ g)„p”°“‘'(^ ; oi, 02 , 03 , 04 !^). 


These orthogonal polynomials are deformation of the Hermite, Laguerre and Jacobi 
polynomials, respectively. 

3.3. Eigenfunctions of single-particle quantum mechanics 

The Hamiltonian m for the single-particle case (re = 1) reads 


H = 


me" 


■^\/wix) e ^medx wix)* a/ wix)* a/ wix) — wix) — wix)*'^ , (34) 


where x = qi. By introducing a dimensionless variable y (or z), a rescaled potential 
Viy) (or Viz)), and a set of parameters A (and a parameter q) given by 


(i) 


A = (01,02), !/(?/) = H(?/; A) = (oi-Fzy)(02-hiy), (35a) 


me 

V = -^x . 

^ h 

TTIC 

y = —x, A = (01,02,03,04), 


Viy) = Viy ^X) = 


111 


Viz) = Viz-,X,q) = 


(oi -F iy)ia2 + iy)ia 3 + iy)iai + iy) 

2iyi2iy + 1) 

A = (oi, 02, o 

(1 - ai2;)(l - a2z)il - a32;)(l - a^z) 


(35b) 


y = —x, z = e y = e l, A = (oi, 02, 03, 04), q = e rncL ^ 


(1 — z‘^)il — qz^) 


(35c) 
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the function w{x) takes the form 




dld2 


huJi ’ hui2 '' ’ 


(ii) : w{x) = -^V{y,\), A = 51,52 + , 

0102 V moi ni 02 / 


huJi ’ 1tJjJ2 


(hi) : i(;(x)* = (010203045 ^^V{z]X,q), X = ,-q^'\. 

Then the Hamiltonian (EH) reduces to the dimensionless one 7i, with 


TTIC 

(i) : H = - Ti , where Ti is given by (|A-18I) , 


11 


H = 


0102 

4mc^ 

01O2 


7i , where 7i is given by (IA-181) , 


(hi) : H = 7710^(010203045 ^) •^Tt, where TL is given by (IA-22j) . 


(36a) 

(36b) 

(36c) 

(37a) 

(37b) 

(37c) 


These Hamiltonians have two properties, factorization and shape invariance 
In the notation used in the appendix, the necessary data for the shape invariance 
are the ground state wavefunction (pQ, the elementary shift 5 of the parameters and 
the first energy level ,?i(A): 


(i) : 

(ii) : 


(hi) : 


(/)o(5; A) oc |r(ai+zy)r(o2+ 75)! , d = (i,i), .Si(A) = Oi + 02 , ( 38 a) 

r(oi+7y)r(o2 + 77/)r(a3 + 7y)r(o4 + 75) 

0 o(y; A) oc 


r{2iy) 

^ = (^, h h l)> + 02 + 03 + 04) , 

A t \ \ iz‘^',Q)c 

(Po{z-,X,q) oc 


(oi 2 ;, 022 :, 032 ;, 042 :; 5 )c 

(5 = ^, (5' = -^, ,5i(A, 5) = i(5"^ - 1)(1 - 01020304). 


(38b) 


(38c) 


From the formulas (lA- 16|) -( E- 17|) / (IA-3 7l)- (|A-38l) and (|A-11|) /(| A-32I) . we obtain the 
eigenfunctions Pn(po and the corresponding eigenvalues Sn of 7i: 


(i) : 

Pniy;X) oc Pniy]ai, 02 , 01 , 02 ), SniX) = ln{n + 2 oi -h 2 o 2 - 1), 

(39a) 

(ii) : 

Pn{y,>^) oc iy„(7/^; 01,02,03,04), 



£n{X) = \n{n -1- oi -1- 02 -1- 03 -1- 04 - 1), 

(39b) 

(hi) : 

Pn{z-, A, 5) oc Pn{P(iZ ] Oi, O2, O3, O4I5), 



£n{X, 5) = \{q~'^ - 1)(1 - 010203045”"^), 

(39c) 


where the first p„, and the second are the continuous Hahn, Wilson and 
Askey-Wilson polynomials in the notation used in Ref. 11). 

We thus find that the equilibrium positions of the multi-particle classical systems 
and the eigenfunctions of the single-particle quantum systems are described by the 
same orthogonal polynomials, namely, the continuous Hahn (special case), Wilson 
and Askey-Wilson polynomials. 
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§4. Summary and comments 


We have established an interesting property of the Ruijsenaars-Schneider-van 
Diejen systems with the rational/trigonometric potentials associated with the clas¬ 
sical root systems. The equilibrium positions of the multi-particle classical systems 
and the eigenfunctions of the single-particle quantum systems are described by the 
same orthogonal polynomials; the continuous Hahn (special case), Wilson and Askey- 
Wilson polynomials. This property is inherited from the Calogero-Sutherland-Moser 
systems, in which the relevant orthogonal polynomials are the Hermite, Laguerre 
and Jacobi polynomials. These polynomials are members of the Askey-scheme of 
the basic hypergeometric orthogonal polynomials. The CSM and RSvD systems ad¬ 
mit elliptic potentials and finding eigenfunctions of such elliptic systems is a good 
challenge. If this property is inherited by the elliptic RSvD systems, study of classical 
equilibrium positions may shed light on the quantum problem of finding eigenfunc¬ 
tions, which is quite non-trivial. 

This interesting property is obtained as a result of explicit computation and at 
present we do not know any deeper reason or meaning behind it. In §6 of Ref. 9), we 
presented the ‘phenomenological’ observation that the Hermite, Laguerre and Jacobi 
polynomials of degree n can be obtained from the prepotential of CSM systems 
with n -|- 1 particles. Is this also true for the continuous Hahn, Wilson and Askey- 
Wilson polynomials? Although there is no prepotential for the RSvD systems, the 
corresponding quantity may be logi^o- The ground state eigenfunction (poiiQj}) of 
(see (Uni) is 


(i) : (poiiqj} ] X, g) oc 


(ii) : 4>o{{qj} ; X, go) (X 


Y[r{ai+iyj)r{a2 + iyj) ■ 


j = l 

n y-fA 

j 

i=i 


n n 


r{g + i{yj -yk)) 
riiigjj - yk)) 

(40a; 

r{go + i{yj + eyfc)) 


(hi) : (l)o{{qj}-,X,go,q) (X 


r{2iyj) 

{zf,q)o, 


l<j'<fc<n €=±1 


n n 


r{i{yj + eyk)) 


izjzl-,q)c 


H Ui=i{aaZj;q)oc i(^0Zjzl-,q) 


(40b) 
, (40c) 


where A is given in (IJ5al) - (IJ0cl) and yj = ^qj-, Zj = and oq = For 

the CSM systems, the second factor of (^o is simply qj ± qj or sm{qj ± qj). The 
observation of Ref. 9) is that the first factor of 4>q gives the measure and the second 
factor gives the desired polynomial. For the RSvD systems, the first factor of (po gives 
the measure, but the second factor does not give the desired polynomial. Therefore, 
the observation of Ref. 9) does not apply to the RSvD systems in its naivest form. 

In the treatment of the Hamiltonians of these single-particle quantum mechanics, 
we have emphasized factorization, shape invariance and construction of the isospec- 
tral Hamiltonians. Although the examples given in this article possess rational and 
trigonometric potentials, this method and the basic idea can be applied to a wider 
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class of potentials, e.g. elliptic potentials. In ordinary quantum mechanics there ex¬ 
ists Crum’s theorem,which describes the construction of the associated isospec- 
tral Hamiltonians Tig and their eigenfunctions 4>s^n for general systems without shape 
invariance. The construction of Tig and 4>g^n given in this article for the ‘discrete’ 
cases is based on shape invariance. A ‘discrete’ analogue of Crum’s theorem, that is, 
the construction of the associated isospectral Hamiltonians and their eigenfunctions 
without shape invariance, would be very helpful, if it exists. 
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Appendix 

Shape Invariance in Single-Particle (‘Discrete’) Quantum Mechanics 


In this appendix we collect useful techniques for the shape invariant single¬ 
particle (‘discrete’) quantum mechanics along the idea of Crum.^®) All the quantities 
in this appendix are dimensionless. 

A.l. Ordinary quantum mechanics 

Let us consider an ordinary quantum mechanical system with the single coordi¬ 
nate y governed by the Hamiltonian 

where the potential V may contain a set of parameters A. We assume that 7i has a 
square integrable ground state and a discrete spectrum, 


Hc^n = £n4>n (u = 0, 1, 2, . . .), 0 = Tq < <?! < ‘?2 < ' ' ' , (A-2) 


where the constant term of V is chosen such that To = 0. Since the ground state 
eigenfunction (j)Q has no node, it can be expressed as </>o(?/) = and the excited 

state eigenfunctions are expressed as (pniu) = Pn{y)<t>o{y)- Then the Hamiltonian li 
can be factorized as 


K = .4fe)UO) = i(-^ + 


( dW{y) \^ d?W{y) \ 

\ dy ) dy‘^ / ’ 


A{y) 


1 / . d . dW{y) \ 
\/2V dy * dy ) 


-4(y)t 


1 / . d 
\/2 V ^ dy 


. dW{y) ^ 

dy J 


(A-3) 

(A-4) 


The ground state cpo is annihilated by A, A(po = 0. 

Let us write the parameter dependence explicitly, as A{y) = A{y ; A). We assume 
the following shape invariance 


A{y ; A)A(y; A)^ = A{y ; A -h d)^A(y; A -F d) -F Ti(A), 


(A-5) 
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where S stands for a set of constants indicating the elementary shift in the param¬ 
eters. Isospectral Hamiltonians can be constructed in the following way.*^ Starting 


from .Ao = A, Ho = H and cpo^n = 4>n-, let us define As, Hg and 4>s,n (n > s > 0) 
recursively: 

As+i{y ] X)'^= As{y; X + S), (A-6) 

Hs+iiy ; A) Agiy ; X)Agiy ; A)t + f,(A), (A-7) 

4>s+i,niy ; A) Agiy ; X)(j)g^niy ; A). (A-8) 

As a consequence of the shape invariance (EiSl), we obtain for n > s > 0 

Agiy ] X) = A{y; X + sS), (A-9) 

; A) = Agiy ; A)^As(y; A) SgiX) = Hiy ; A s5) + ^^(A), (A-10) 

n—1 

Sg+iiX)=SgiX) + SiiX + sS) [=^SniX) = ^SiiX + sS)), (A-11) 

s=0 

T~(-siy ; X)(j)g^niy ; A) = SniX)(j)g^niy ; A), (A'12) 

Agiy; X)(j)g^giy ■ X) = 0, (A-13) 

Agiy ; A)Vs-hi,n(y; A) = (<?n(A) - <5*(A))(y ; A). (A-14) 

As is clear from (lA-lOll . all the Hamiltonians Hq, Hi, ..., Hg, ... have the same shape 


and only the parameters are shifted. From (EiHl) and (IXdlll we obtain formulas 
relating the eigenfunctions along the horizontal line (the isospectral line) of Fig. 1, 


4>s,niy; A) = Ag-iiy ; A) • • • Aiiy ; A)Ao(y; X)4>niy ; A), 

, . A(y;A)^ Ai(y;A)t A„_i(y;A)t 

^ Sr,iX)-SoiX) SniX)-SiiX)'" 
and from (jA-IDI) we have 


£niX) - £n-liX) 


4^n,niy : A) , 


(A-15) 

(A-16) 


(l)n,niy;X) (X (l)oiy;X + n6). (A'17) 

It should be emphasized that all the operators A and A^ in the above formulas are 
explicitly known thanks to the shape invariance. The formula (jA-l611 with (IA-1 7|1 
can be understood as the Rodrigues-type formula. The relation (lATlI) implies that 
{T„(A)}n>o can be calculated from Ti(A), and hence that the spectrum is determined 
completely by the shape invariance. 

As seen above, the operators A and A^ act isospectrally, that is, horizontally in 
Fig. 1. On the other hand, the annihilation and creation operators map from one 
eigenstate to another (i.e. vertically) of a given Hamiltonian. In order to define the 
annihilation and creation operators, we need unitary operators lAg. (For details see 
for example Refs. 18),20),28).) 

The scheme described above is illustrated in Fig. 1. 

Isospectral Hamiltonians can be constructed without shape invariance and some formulas e.g., 
and lEdeJ, are also valid in that case. However, there is no explicit systematic method to 
obtain .As (and £n) in general. 
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^3 


£2 






Fig. 1. A schematic diagram of the energy levels and the associated Hamiltonian systems together 
with the dehnition of the A and operators and the ‘creation’ (a^) and ‘annihilation’ (a) 
operators. The parameter set is indicated below each Hamiltonian. 

A.2. ‘Discrete’ quantum mechanies 

Here we discuss two different types of shape invariance mechanisms. The first, 
with an additive shift of the parameters, appears in the ordinary quantum mechanics 
discussed in the preceding subsection. The second one, involving a multiplicative 
shift of the parameters, is new. 

A.2.1. Additive shift 

First, let us consider the Hamiltonian 

^ \ (y) (y)* + (y)* ^y) ~ ^(y) ( ats ) 

and its eigenfunctions given in (IX^ . The function V(y) depends on a set of param¬ 
eters A. This Hamiltonian TC is factorizable; 

n = n{y,X)=A{y,X)U{y,X), (AT9) 

A = A{y ; A) ^ {y ; A)* - ^ ^JV (y ; A) ^ , (A-20) 

A^ = A(y; A)t ='" ■^(^^/V{y,X) e~^^ - y/V{y]X)* e^^'j . (A-21) 

The ground state 4>o is annihilated by A, Acpo = 0. 

As in go we assume the shape invariance (IOi; that is, AA"^ has the same 
shape as A'^A with the additive shift of the parameters A combined with the ex- 
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citation energy of the first level fi(A). Then the isospectral Hamiltonians can be 
constructed in the same way as in the previous subsection and equations (|A-6|) 
also hold here.^^^ 

A.2.2. Multiplicative shift 

Next let us consider the Hamiltonian 

n = ^(^^/V{^)q^^^/vW+VWrQ~'^^Vnz)-V{z)-V{zY) (A-22) 

and its eigenfunctions (IX^ . The variable z is related to y as z = and = z^. 
(In our notation, V{z) here corresponds to V{y)* in (IA-18|1 .1 This form is more 
convenient for the trigonometric potentials. The function V{z) depends on a set 
of parameters A and a parameter q, which controls the multiplicative shift. This 


Hamiltonian TC is factorizable: 

Ti. = TC{z ]X,q) = A{z ; A, q)'^A{z ]X,q), (A-23) 

A = A(z ; A, q) = ^ [q-2^^ ^V{z-X,qY - ) , (A-24) 

= Aiz-,X,qY Viz] X,q)q^^^ - y/Viz ] X, q)* . (A-25) 

The ground state cpo is annihilated by A, Acpo = 0 . 

We assume the following type of shape invariance 

Aiz ; A, q)Aiz ; A, = q^^'Aiz ; q^X, qYAiz ; q^X, q) + £iiX, q ), (A-26) 

where 6 and 6' are constants. As in go isospectral Hamiltonians can be con¬ 
structed.^^) Starting from Aq = A, Hq = H and (^o,n = (t>n, let us define Hs 
and (j)s^n (n > s > 0) recursively as 

As+iiz ; A, q) q^'Asiz ; q^X, q ), (A-27) 

Hs+iiz ; A, q) Asiz ; A, q)Asiz ; A, qY + ^^(A, q ), (A-28) 

4>s+i,niz ; A, q) Asiz ; A, q)(j)s,niz ]X,q). (A-29) 

As a consequence of the multiplicative shape invariance (|A-26I) . we obtain for n > 
s > 0 

Asiz ]X,q) = q'^^'Aiz ; q^^X, q), (A-30) 

TCsiz ]X,q) = Asiz ; A, qYAsiz ]X,q)+ Ts(A, q) 

= q^'^^'niz ; q^^X, q) + Ts(A, q ), (A-31) 

n—1 

£s+iiX, q) = £siX, q) + q^^^'S^q^'X, q) T„(A, q) = Y. q^^^'S^q^'X, q)) , (A-32) 

s=0 

Hsiz ] A, q)(j)s,niz ]X,q) = T„(A, q)(j)s,niz ]X,q), (A-33) 

Asiz]X,q)(j)s,siz]X,q) = 0, (A-34) 

Asiz ] A, qY(j)s+i,niz ]X,q) = (Tn(A, q) - Ts(A, q))(j)s,niz ]X,q). (A-35) 
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It should be noted that the Hamiltonian 7i is rescaled by a factor q'^^' and the 
parameters A are multiplied by a factor at each step. From (IA-291) and (|A-35ll we 
obtain the formulas 


4>s,n{^ ■, X, q) ■ ■ ■ Ai{z ] q)Ao{z ; A, q)(l)n{z ;X,q), 

Ao{z;X,qy Aiiz;X,q)^ 


(A-36) 


i{z',X,q) = 


Sn{X, q) - So{X, q) £’„(A, q) - Si{X, q) 
^ _ An-i{z]X,qy 

£n{X,q) - £n-i{X,q) 


4>n,n{^X, q) , (A-37) 


and from (IA-311) we have 


4>n,n{^ ; 9) oc (l)o{z ] q'^^ X, q). 


(A-38) 
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